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ABSTRACT 

We present a population study of several types of galaxies within the protocluster surrounding 
the radio galaxy MRC 0316-257 at z ~ 3.1. In addition to the known population of Lya 
emitters (LAEs) and [Oin] emitters, we use colour selection techniques to identify protocluster 
candidates that are Lyman break galaxies (LBG) and Balmer break galaxies (BBGs). The 
radio galaxy field contains an excess of LBG candidates, with a surface density 1.6±0.3 times 
larger than found for comparable blank fields. This surface overdensity corresponds to an 
LBG volume overdensity of ~ 8 ± 4. The BBG photometric redshift distribution peaks at the 
protocluster' s redshift, but we detect no significant surface overdensity of BBG. This is not 
surprising because a volume overdensity similar to the LBGs would have resulted in a surface 
density of ~ 1 .2 that found in the blank field. This could not have been detected in our sample. 
Masses and star formation rates of the candidate protocluster galaxies are determined using 
SED fitting. These properties are not significantly different from those of field galaxies. The 
galaxies with the highest masses and star formation rates are located near the radio galaxy, 
indicating that the protocluster environment influences galaxy evolution at z ~ 3. We conclude 
that the protocluster around MRC 0316-257 is still in the early stages of formation. 

Key words: galaxies: evolution - galaxies: high-redshift - galaxies: clusters: general - cos- 
mology: observations - cosmology: early Universe 



1 INTRODUCTION 

One of the main aims of astrophysics is to understand the formation 
and evolution of galaxies. Hierarchical evolution in ACDM cos- 
mology means the evolution of a galaxy will depend on whether 
it is l ocated in a low-density or high-density environment i Toomrd 
Il977|). This has been quantified in studies bv lClemens et alJfeOOfj) 
and Sanchez-Blazque z et al.l d2006h for the local Universe and for 
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the more distant Univer se bv lvan Dokkum & van der Marell ( 2007) 
and iGobat et"al] J2008h . These studies find that early-type galax- 
ies in cluster environments are older than early-type galaxies resid- 
ing in low-density, field environments, indicating that galaxies in 
cluster-like environments form at an earlier epoch. 

Further evidence that environment influences galaxy evolu- 
tion is the observation of ' environmen t-dependent downs izing ' in 
galaxy clusters at z ~ 1. Downsizing dCowie etalJll996l) implies 
that the bright, massive galaxies move onto th e red sequence first , 
whilst fainter galaxies are added at a later time. lTanaka et ail feooi 



2 E. Kuiper et al. 



I200I I2008I) show that the red sequence in high density environ- 
ments extends to fainter magnitudes than in less dense galaxy 
groups. The red sequence therefore forms at an earlier epoch in 
dense environments. 

To understand galaxy evolution in different environments it 
is necessary to study galaxy clusters across cosmic time. Galaxy 
clusters have been detected out to z = 1.5 using conventional tech- 
niques such as through observation of the hot X -ray emitting intra- 
cluster gas or IR red sequ ence searches (e.g.. iMullis et al. 2005; 
IStanford let alj2005ll2006h . 

The most successful technique for finding cluster progenitors 
at z > 1.5 is to search for emission line galaxies around high- 
redshift radio galaxies (H zRGs, for a comprehensive review see 
iMilev & De Breuckll2008h . HzRGs are among the most luminous 
and massive objects in the early Unive rse and are expected to 
be the progenitors to local cD galaxies dRocca-V olmerang e et al.l 
|2004| ; ISevmour et alj[2007h . Several studies have found that they 
are situated in overd ense regions with prop e rties expected of form- 
ing galaxy clusters dPentericci et all [2 000 : Venemans et alj[2005L 



120071 ; llntema et alj|2006l ; bverzier et alj|2006l . 120081) . The dense 
cluster-like environments around HzRGs are likely not yet virial- 
ized and are commonly termed 'protoclusters'. They are excellent 
laboratories for studying the formation and evolution of galaxies in 
overdense environments. 

In this study we investigate galaxy populations around the 
HzRGMRC 0316-257 (hereafter 0316) located atz = 3.13. Galaxy 
selection techniqu es based on the Lyman break at 912 A (e.g., 
ISteidel et alj|2003l) are most efficient at this redshift. Also at this 
redshift, strong emission lines fall within existing narrowband fil- 
ters. This allows the selection of many different types of emis- 
sion line galaxies in the protocluster. Therefore the redshift of 
this protocluster make s it ideal for studying its galaxy populations. 
Venemans et al. I d2005h (hereafter V05) sh owed this region contains 
an overdensity of hya emitters (LAEs) and|Maschietto et alj J2008h 
(hereafter M08) found a number of [Om] emitters near the redshift 
of the radio galaxy (RG). Additional galaxy populations are iden- 
tified using a large set of broadband images. The properties of the 
galaxy populatons are determined using broadband photometry and 
spectral energy distribution (SED) fitting. These properties are then 
compared to the properties of field galaxies to search for environ- 
mental influences on galaxy evolution at this redshift. 

The paper is ordered as follows: in Sect.|2]the data and its re- 
duction is discussed. This is followed by the photometry in Sect. [3] 
and sample selection in Sect. [4] The particulars of the SED fitting 
process are discussed in Sect. [5] after which we show our results 
in Sect. [6] We discuss and compare the results to work by other 
authors in Sect. [7] Finally, the summary and conclusions are pre- 
sented in Sect. [8] Throughout this paper we use a standard ACDM 
cosmology with H =71 km s" 1 Mpc" 1 , fi M =0.27 and fi A =0.73. All 
magnitudes given i n this paper are in the AB magnitude system 
(Oke &Gunil 19831) unless noted otherwise. 



2 DATA 

Images of the 0316 field were obtained in 20 passbands spanning 
the U band to the 8 fim band. A summary of all the data used is 
given in Table 1 and the filter response curves of each of the filters 
is shown in Fig. [T] togethe r with two example SEDs taken from the 
iBruzual & Charlotl fe003h (BC03) models. The fields covered by 
the various instruments are illustrated in Fig. [9] The reduction of 
the various data sets is described below. 



2.1 Ground-based UV-optical imaging 

UBVR imaging data were obtained using the Visible MultiOb- 
ject Spectrograph (VIMOS. lLe Fevre et alj|2003ri instrument at the 
Very Large Telescope (VLT). The data were taken during the pe- 
riod of 14-15 November 2003 for the V band and 20-25 Novem- 
ber 2003 for the remaining bands. Additional U band data were 
obtained on 14 February 2004. 

The VIMOS field-of-view consists of 4 separate quadrants, 
each having a field-of-view of approximately 7'x8'. Two pointings 
were used in which the RG was centered in one of the quadrants. 
For most of this work only the central quadrant that contains the 
RG is used, as only this central region has additional data. 

The reduction of the UV and optical data was performed us- 
ing standard tasks in the IRAff^l software package. The process 
includes bias subtraction and flat fielding using twilight sky flats. 
Remaining large scale gradients were removed using a smoothed 
master flat. This was obtained by median-combining the unregis- 
tered, flatfielded science images. Reduced images were registered 
and combined to form the final science images. Photometric ze- 
ropoints were determined using standard star images taken on the 
same nights as the science frames. 

Deep Keck u' band data (14) were obtained on 31 January, 1 
February and 4 February 2003 using t he blue arm of t he Low Reso- 
lution Imaging Spectrometer (LRIS. IOke et aljfl995h . As there are 
large differences in the filter responses, both U bands have been in- 
cluded in the analysi s. Information regardin g the reduction of these 
data can be found in lVenemans et alj d2007h . 

V and / band data of the 0316 field were obtained during the 
respective periods of 20-21 and 6-8 September 2001 using the 
FORS2 instrument at VLT. Information on these data and their re- 
duction can be found in V05. 



2.2 HST/ACS optical imaging 

Deep Hubble Space Telescope (HST) ^25 and Igu images of 
the 0316 field cover approximately half of the VIMOS field-of- 
view. These images were obtained using the Advanced Camera 
for Surveys (ACS, iFord et alj Il998h during the periods of 14- 
31 December 2004 and 2-21 January 2005. The data covered 
two fields of 3.4'x3.4' with approximately 1' overlap between the 
two fields. The two / 814 fields were combined with an additional 
3.4'x3.4' ACS field which was obtained on 18 July 2002. Details 
of the data and their reduction can be found in M08 and V05. 



2.3 Near-infrared data 

Two sets of near-infrared (NIR) data spanning the J to K s bands 
were used. Deep J and A" s images were ob tained with the Infrared 
Spectrometer And Array Camera (ISAAC, iMoorwoodelafll 19981) 
on the VLT on various dates between November 2003 and October 
2004 (see M08). These images are deep, but only cover the inner- 
most 2.5' x2. 5' of the protocluster. 

Additional JHK^ data were obtained using th e Multi- Object 
InfraRed Camera and S pectrograph (MOIRCS, llchikawa et al.1 
2006; Suzuki et alj|2008h at the Subaru telescope on 6-7 January 
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omy, Inc., under cooperative agreement with the National Science Founda- 
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Table 1. Details of the observations. The 5cr limiting magnitudes have been calculated for angular diameters of 2" for the ground based data, 0.5" for the ACS 
data and 4" for the IRAC data. 
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Instrument 


Date ol observation 


A eff( A > 




Exp. time (sec.) 


Seeing 


5o~ limiting mag. 




LRIS/Keck 


2003 Jan. 31, 2003 Feb. 1 &4 


3516 


561 


17100 


1.3" 


26.33 




VIMOS/VLT 


2003 Nov. 23-24, 2004 Feb. 14 


3744 


359 


15500 


1.0" 


25.98 


B 


VIMOS/VLT 


2003 Nov. 22 


4310 


832 


3720 


0.8" 


26.02 




VIMOS/VLT 


2003 Nov. 14-15 


5448 


842 


4960 


0.75" 


25.97 


v t 


FORS2/VLT 


2001 Sept. 20-21 


5542 


1106 


4860 


0.7" 


26.10 


R 


VIMOS/VLT 


2003 Nov. 20 & 25 


6448 


1292 


6845 


0.9" 


25.90 


1 


FORS2/VLT 


2001 Sept. 6-8 


7966 


1433 


4680 


0.65" 


25.83 


r 625 


ACS/HST 


2004 Dec. 14-31, 2005 Jan. 2-21 


6321 


1327 


23010 


- 


27.66 


h u 


ACS/HST 


2002 Jul. 18, 2004 Dec. 14-31, 2005 Jan. 2-21 


8089 


1765 


52320 




28.37 


h 


ISAAC/VLT 


2003 Nov. - 2004 Oct. 


12535 


2640 


19000 


0.5" 


24.60 




ISAAC/VLT 


2003 Nov.- 2004 Oct. 


21612 


2735 


19000 


0.5" 


24.19 


J m 


MOIRCS/Subaru 


2006 Jan. 6-7 


12532 


1538 


4680 


0.75" 


23.77 


H 


MOIRCS/Subara 


2006 Jan. 6-7 


16364 


2788 


3600 


0.8" 


22.71 


Urn 


MOIRCS/Subaru 


2006 Jan. 6-7 


21453 


3042 


3300 


0.75" 


23.07 


[3.6] 


IRAC/Spitzer 


2005 Jan. 17-19 


35636 


6852 


46000 




23.76 


[4.5] 


IRAC/Spitzer 


2005 Jan. 17-19 


45111 


8710 


46000 




23.60 


[5.8] 


IRAC/Spitzer 


2005 Jan. 17-19 


57598 


12457 


46000 




22.24 


[8.0] 


IRAC/Spitzer 


2005 Jan. 17-19 


79594 


25647 


46000 




21.97 


Lya NB 


FORS2/VLT 


2001 Sept. 20-21 


5040 


61 


23400 


0.7" 


25.3 


[Ora] NB 


ISAAC/VLT 


2003 Nov. - 2004 Oct. 


20675 


437 


24840 


0.45" 


22.6 



2006. For d etails concerning the reduction of these data we refer 
the reader to lKodama et all d2007ri . hereafter K07. This set of JHK S 
imaging data is shallower than the ISAAC data described above, 
but it covers a larger fraction of the 0316 field. To avoid confusion 
between the ISAAC and MOIRCS data we denote the bands with a 
subscript T or 'm', respectively. 

2.4 Mid-infrared data 

In the mid-infrared w avelength range w e have Spitzer InfraRed 
Array Camera (IRAC. iFazio et alj|2004h data at 3.6 yum, 4.5 /xm, 
5.8 fim and 8.0 yum (hereafter [3.6], [4.5], [5.8] and [8.0] respec- 
tively). IRAC data were obtained in all bands on 17-19 January 
2005 covering a ~5'x5' field centred on the RG. Deep imaging 
was obtained using a medium scale cycling dither pattern of 230 
frames with a 200 s frame time for a total exposure time of ~12.7 
hours. The [3.6] and [4.5] basic calibrated leve l data (BCD) were 
reduc ed and mosaiced using the mopex software ( Makovoz & Khan! 
l2005h following standard procedures. Before the BCD frames were 
combined, the muxbleed and column pulldown effects were cor- 
rected using custom software provided by D. Stern and L. Mous- 
takas. During mosaicing, the images were resampled by a factor 
of V2 and rotated by 45°. The [5.8] and [8.0] BCD data were fur- 
ther mosaiced using a custom IDL code kindly provided to us by I. 
Labbe. 

2.5 Further reduction 

With the exception of the ACS data, all images were resampled to 
a common pixel scale of 0.205"/pixel and transformed to the same 
image coordinate system using the IRAF tasks geomap and greg- 
ister. The images were then convolved with 2D Gaussian profiles 
to match the PSF FWHM of the VIMOS U band (l/ v ), which has 



the largest seeing of approximately 1". The IRAC bands and the 
{7 k band are excluded from this process as they have significantly 
larger PSF sizes. Smoothing the other images to the PSF size of 
these images would negatively impact the quality of the analysis. 

Because of the extreme smoothing required to match the res- 
olution of the ACS images to that of the ground-based data, the 
ACS data were not used for the SED fitting and photometric red- 
shift determination, but only for the determination of UV slopes of 
the protocluster galaxy candidates. 



3 PHOTOMETRY 

Photometry was obtai ned using the SExtractor software 
dBertin& Amoutsl Il996h in double image mode. Lyman Break 
Galaxy candidates (LBGs) were detected using the unsmoothed R 
band as the detection image. For the Balmer Break Galaxy can- 
didates (BBGs) the unsmoothed K s bands were used as detection 
images. 

A detection is defined as an object with 5 adjacent pixels that 
each exceed the 3cr rms noise. Colours were measured in the 3cr 
level isophot apertures as determined from the detection image. If 
a known object was undetected in an image it was assigned the 3cr 
detection limit of the band in question. Lowering this limit to lcr 
does not affect any of the conclusions presented in this paper. To- 
tal magnitudes and fluxes in the detection bands were obtained by 
using the MAG_AUTO function of SExtractor. Total fluxes in the 
remaining bands were obtained by scaling the respective isophotal 
magnitudes accordingly. 

Because of the large size (FWHM) of the PSF in the IRAC im- 
ages it was not possible to determine consistent colours using the 
isophotal apertures mentioned above. Instead, the object flux was 
measured in a circular aperture of 20 pixel (4"), and a correction 
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3000 4 000 5000 6000 7000 8000 9 000 2x10 4 4x10" 6x10 4 8x10" 1x10 5 

Wovelength(Angstrom) Wovelength (Angstrom) 



Figure 1. The filter response curves for the filters used in this study. For illustrative purposes the curves have been scaled to the same maximum throughput. 
The J m and K m filter curves are not shown due to their similarity to the corresponding ISAAC filters. Also plotted are two z = 3.13 model SEDs obtained from 
the BC03 population synthesis models. The light gray SED is for a continuous star forming galaxy of 100 Myr old, whereas the dark gray SED is for a 1 Gyr 
old galaxy which has an exponentially declining star formation with r = 10 Myr. 



factor applied to match the photometry to the other bands. This cor- 
rection factor was determined by smoothing the detection images 
to the spatial resolution of the IRAC data. Object fluxes were mea- 
sured in the smoothed detection images in 4" circular apertures, 
and compared to the flux measured from the unsmoothed detection 
images in the standard 3cr isophotal apertures. The ratio between 
the two fluxes yielded the correction factor. A similar process was 
applied to determine the colours in the Uy, image. 

The large FWHM of the PSF in the IRAC bands also causes 
source confusion and contamination by neighbouring sources. For 
faint objects this effect can be a large source of error. All detected 
objects were visually inspected for contamination, and heavily con- 
taminated objects have been removed from the analysis when rele- 
vant. 

Finally, all magnitudes were corrected for Galactic foreground 
extinction determined from the Schlegel et al. (1998) extinction 
maps. 

Photometric uncertainties and limiting magnitud es (listed in 
Table TJ were computed using the method described in Labbe et al. 

which is summarized below. After masking all objects, the 
rms of pixels in the entire image and the rms of fluxes in aper- 
tures of various sizes were measured. We determined the noise in 
an aperture of size N using the relation 



o-,(N) = N&iidi + hN) 



(1) 



with <fi the rms of pixels over the entire image and cr, the rms within 
a certain aperture size N. N is defined as N = VA where A is 
the area of the aperture. The subscript T indicates the photomet- 
ric band in question. The free parameters a.\ and b, were then fitted 
such that it can be calculated what the noise is in an aperture of a 
given size. The uncertainty calculated using the rms of all the back- 
ground pixels can be an underestimate, because it does not take into 
account pixel-to-pixel dependencies introduced in the reduction of 
the data. Due to the IRAC photometry being more uncertain an ad- 
ditional 10 per cent uncertainty was added in quadrat ure for those 
four b ands as has been done in previous studies (e.g. lLabbe et al.l 
l2005t) . 

The completeness of the detection fields for point sources is 
shown in Fig. [2] The completeness was determined by extracting 
several bright, unsaturated stars from the images, averaging them to 
obtain a high signal-to-noise PSF image and then inserting a num- 
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Figure 2. Completeness of the three detection bands as a function of mag- 
nitude for points sources. 



ber of these PSFs at random locations in the images. The source 
extraction was repeated and the number of recovered artificial stars 
yielded the completeness as a function of magnitude. To avoid over- 
crowding only 150 objects were added at a time. This process was 
repeated ten times for each magnitude to obtain better statistics. 
Due to the small field size of the ISAAC data the number of added 
objects was lowered to 50 and the process was repeated fifty times. 
The data are 50 per cent complete down to R = 26, K\ = 24.1 and 
K m = 23.2. 



4 SAMPLE SELECTION 

4.1 Lya- and [Om] excess objects 

Venemans et al. I J20051) spectroscopically confirmed 32 LAEs at the 
redshift of the RG, and found that the volume density of LAEs in 
the 0316 field is 2-4 t i mes la rger than for the blank field. 

Maschietto et al. (2008j) found a sample of 13 [Om] emitting 
galaxies near the RG, corresponding to a surface density 3.5^t 
times the field density (uncertainties obtained from Poisson statis- 
tics). The [Om] overdensity is consistent with the LAE overdensity 
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found by V05, but the sample is small and the uncertainties are 
large. Five of the [Om] emitters are spectroscopically confirmed 
LAEs, and three additional [Om] emitters are spectroscopically 
confirmed to be at z = 3.1. These [O m] emitters are blueshifted 
with respect to the RG indicating that the protocluster around 03 16 
may be part of a larger superstructure. 



4.2 Lyman Break candidates 

LBGs were selected using a colour criterion similar to that used by 
ISteidel et alj d2003h to select star forming galaxies at z ~ 3. The 
criterion uses the U^VR bands. Even though the VIMOS field is 
larger, the U y passband is redder than if passbands used in most 
other studies. Thus the 14 band was used to facilitate comparison 
with other LBG studies. 

The colour criterion was devised by creating artificial galaxy 
spectra using the BC03 evolutionary population synthesis models. 
Galaxy spectra were synthesized for various star formation histo- 
ries (SFHs) and a variety of values for extinction and age. The spec- 
tra were redshifted from z = to z = 5. The model spectra were 
then convolved with the filter curves to obtain synthetic galaxy pho- 
tometry. Galaxies situated at 3.0 < z < 3.3 lie in the upper left cor- 
ner of the £4 - V v vs. V v - R colour-colour diagram. This region is 
parametrized by the relations 



t/ k - y v > 1.9, 

V v -R<0.5\, 
U k -V v > 5.07 x (V v -R) + 2.43, 
R< 26. 



(2) 



The £4 - V v vs. V v - R colours for the R band detected sam- 
ple are shown in Fig. [3] The LBG selection criterion is marked as 
solid lines. The inset shows the ratio of the number of synthesized 
objects that are selected to the total number of synthesized objects 
as a function of redshift. This criterion should predominantly se- 
lect objects having redshifts between 2.9 < z < 3.4 which have low 
ages and little dust obscuration (t < 100 Myr, E(B - V) < 0.3). The 
50 per cent completeness limit of R = 26 was adopted as a magni- 
tude cut. A total of 52 objects in the 0316 field satisfy the selection 
criterion. 

Photometric redshifts (z p hot) were det ermined for all R ban d 
detected objects using the EAZY code (Bra mmer et al.1 2008). 
Since a significant number of objects lack deep NIR coverage we 
include the more uncertain IRAC bands. To test the influence of 
the IRAC data on the photometric redshift determination we have 
determined the photometric redshifts both including and excluding 
the IRAC photometry. 

The EAZY code yields 4 different estimates of the redshift. 
It offers the options of applying both marginalizatior0 and a 
Bayesian prior. In order to determine which option is best suited 
for this study we used the Multiwavelength Survey by Yale- 
Chile (MUSYC) ECDF- S data dGawiser et al.ll2006bl : lDamen et al.1 
l2009l ; lTavlor et al]|2009h to produce an R band detected catalogue. 
The photometric redshifts of this sample were determined using 
EAZY and subsequently compared to a lar ge sample of spectro - 
scopic redshifts (z SDec ) fro m the catalogues oflCimatti et alj d2002h . 
iLe Fevreetail d2004l) and lRavikumar et al.1 d2007l) . A total of 2469 
objects have spectroscopic redshifts. Note that this will not be a 



14 The process of marginalization calculates the best photometric redshift 
value by weighting it according to the redshift probability distribution. 




Figure 3. £4 - V* v vs. V v - R diagram of all objects detected in the R band. 
Diamonds indicate the points with lower limits to their £4 _ V v colours. 
The lines indicate the area in colour-colour space that was used for the LBG 
selection. The inset shows the selection efficiency for all objects (full line), 
objects that are younger than 100 Myr and have E(B-V) < 0.3 (dotted line) 
and objects that are older than 500 Myr (dashed line). The criterion should 
predominantly select objects in the range 2.9 < z < 3.4. 
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Figure 4. z p hot versus z spcc for a sample of 2469 objects detected 
in the ECDF-S. No magnitude cut was applied. Photometric redshifts 
were obtained using the EAZY redshift code and include data in the 
UBVRIJHK S [3.6][4.5] bands. The four panels represent the results for the 
four different redshift determinations. Each panel indicates whether a prior 
or marginalization was used. The vertical lines mark the redshift interval of 
interest. It is apparent that it is best to use neither prior nor marginalization 
for z > 2 objects. 



completely consistent test case for the 0316 data, because the ob- 
jects for which spectroscopic redshifts are available are in general 
brighter than R = 25, and the details concerning the depths in cer- 
tain MUSYC bands are different from the 0316 dataset. 

The comparison between z spec and z p h ot for the MUSYC galax- 
ies is shown in Fig. [4] The four panels show the results for four 
different methods of obtaining the photometric redshift. The scat- 
ter is quantified by the median absolute deviation of dz = fcpec _ 
£ P hot)/(l + Zspcc)- For z spcc < 1.5 applying both the prior and the 
marginalization yields the least amount of scatter (\dz\ = 0.043) 
and gives in general good agreement. However, application of the 
prior also shifts the objects with z spec > 2 to z P hot ~ 0.1. 
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Figure 5. Redshift distributions of the sources selected using the colour 
cuts in Eq.|2] The full line and dotted line denote results obtained with and 
without the IRAC bands, respectively. In the inset the redshift uncertainty 
Az is plotted as a function of R band magnitude for all object that satisfy 
the LBG colour selection criterion. Also shown are average values for 0.5 
mag bins. Objects that have uncertainties larger than 1 (as indicated by the 
horizontal line) are not included in these average values. These objects are 
degenerate between z ~ and z ~ 3 and would skew the average to higher 
values. 

Since we target galaxies at z ~ 3, the photometric redshift de- 
termination without prior or marginalization yields the best results. 
Comparing the scatter for z spec > 2 yields \dz\ = 0.22 for the up- 
per left panel and \dz\ = 0.69 and \dz\ = 0.65 for the upper and 
lower right panels, respectively. The option that uses neither prior 
nor marginalization systematically overestimates the redshifts of 
objects at z ~ 1, placing them at 2 < z phot < 2.5 instead. However, 
only 0.2 per cent of objects with z spcc < 1.5 are placed in the range 
2-7 < Zphoi < 3.5, so this should have no significant effect on the 
results. In the remainder of this paper, photometric redshifts of the 
0316 galaxies have been obtained using neither prior nor marginal- 
ization. 

The photometric redshift distribution of the 0316 LBG candi- 
dates is shown in the main panel of Fig. [5] The majority of objects 
is situated at z ~ 3.1 and strong peaks are seen at the redshift of 
the protocluster and z ~ 2.9. These peaks are artefacts caused by 
the discrete number of possible redshifts. Excluding the IRAC data 
(dotted line) has little influence on the photometric redshift distri- 
bution. This is due to the inclusion of both 14 and U v . Both bands 
sample the strong Lyman break feature of these galaxies and hence 
the photometric redshifts are well constrained. 

All objects within 2.7 < z < 3.5 were selected for the fi- 
nal sample of 48 objects. This is based on the photometric redshift 
uncertainties shown in the inset of Fig. [5] Here the redshift uncer- 
tainty Az is taken to be half of the lcr uncertainty interval yielded 
by EAZY. The inset indicates that an uncertainty of 0.4 is a reason- 
able value to encompass all possible z ~ 3.1 objects. The accuracy 
of photometric redshifts is not sufficient to verify if any of these 
objects are indeed members of the protocluster surrounding 0316; 
to achieve this spectroscopic redshifts are needed. 

4.3 Balmer Break Galaxy candidates 

The Balmer break can be used to select z ~ 3 galaxies in a similar 
fashion to the Lyman break method discussed in Sect. 14. 21 

Extracting a sample of distant red galaxies (DRGs) was done 



Figure 6. Redshift distribution of the ISAAC detected BBGs and the 
MOIRCS detected BBGs are shown as a dotted and dashed line, respec- 
tively. The redshift distribution for the combined sample is also shown as 
the full line. Note the peak in number at the redshift of the radio RG. 



by using th e sim ple colour cut (J - /sT 6 )vega S 2.3 proposed by 
iFranx et al] ( 120031) . The DRG colour criterion samples red, pre- 
dominantly massive galaxies between a redshift of 2 and 4. Al- 
though the NIR filter sets used in this study differ from other stud- 
ies the (J - K) colour difference is typically ~0.02 mag, which is 
small compared to the photometric errors. Therefore these colour 
terms are deemed negligible. 

Thirty-four DRGs are found in the ISAAC field-of-view af- 
ter removing 6 ghost images. After applying a 50 per cent com- 
pleteness cut at Ks = 24 the final sample of 17 DRGs remains. 
Amongst the DRGs are the RG, an LAE and an [Om] emitter with 
z spec = 3.104 (M08). 

The MOIRCS data is about 1 magnitude shallower in the K s 
band than the ISAAC data. K07 apply a magnitude cut of K m = 
23.7 and find a total of 54 DRGs in the MOIRCS field, 14 of which 
are located in the ISAAC field. Using the detection criterion de- 
scribed above (5 adjoining pixels 3<r above the background) leads 
to several spurious detections in the MOIRCS data compared to the 
ISAAC data. Increasing the detection criterion to 3.25<r removes 
most of the spurious detections including two DRGs detected by 
K07. These are spurious detections since neither of these two ob- 
jects have counterparts in the deeper ISAAC data. There may also 
be spurious DRGs outside the ISAAC field-of-view. DRGs with 
K m > 23 (the 50 per cent completeness level) were removed from 
the sample resulting in a final sample of 23 DRGs in the MOIRCS 
field. The difference in the number of DRGs in the K07 study and 
the present study is due to the stricter magnitude cut and the differ- 
ent parameter values used for source detection. 

Combining the two separate DRG samples leads to a sample 
of 33 unique DRGs. In the rest of the paper we will refer to DRGs 
as Balmer Break Galaxies (BBGs). 

The main panel of Fig. [6] shows the photometric redshift dis- 
tribution that is found for both BBG samples and the combined 
sample. We find that there is a peak at z ~ 3, the redshift of the 
protocluster. The inset shows the uncertainties on the photometric 
redshifts as a function of K s magnitude. As can be seen from the 
inset the uncertainty on the photometric redshift increases rapidly 
as the magnitude becomes fainter, with lcr uncertainties > 0.5 for 
K\ > 23. This is due to the faint {/ k and U y magnitudes of these 
objects resulting in a poorly constrained Lyman break. Such large 
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Figure 7. z p hot versus z spec for all the objects in the 0316 field that have 
spectroscopic redshifts. The line indicates equality. The main panel shows 
the redshift range of interest, whereas the inset shows low redshift objects. 
Circles denote the 13 objects that are unclassified, the diamonds are LBGs, 
the blue crosses denote the LAEs and the green triangles indicate the [Om] 
emitters. The RG is the [Om] emitter in the lower right corner of the figure. 



Figure 8. Redshift distributions of 'non-classified' objects detected in R 
band that have 2.7 < z p h t < 3.5 and R < 26. The inset shows the subset 
of objectts that satisfy the redshift, V v - R and magnitude criteria for LBG 
selection but are undetected in U^. 



uncertainties would smooth any peak in the redshift distribution 
caused by the protocluster. 

A Kolmogorov-Smirnov test was used to compute the signif- 
icance of t he peak using the ph otometric redshift distribution of 
BBGs from Graz ian et alj d2007l) as a comparison sample. There is 
a probability of 0.3 per cent that both distributions have been drawn 
from the same parent distribution. Thus the photometric redshift 
distribution of BBGs in the 0316 region differs at the 3<r level from 
that in the field. 

Possible protocluster members are identified as the BBGs that 
lie within 2.7 < z phot < 3.5. A total of 9 ISAAC BBGs satisfy this 
criterion. For the MOIRCS BBGs a similar fraction of 1 1 out of 23 
lie within 2.7 < z phot < 3.5. 



4.4 Spectroscopic redshifts 

Within the 0316 field spectroscopic redshifts are known for all 32 
LAEs, 7 [Om] emitters and 16 additional objects (Bram Venemans, 
private communication). Fig. [7J compares these spectroscopic red- 
shifts to the photometric redshifts obtained with the EAZY photo- 
metric redshift code. 

In general the photometric redshifts agree well with the spec- 
troscopic redshifts. Even the intrinsically faint LAEs (denoted by 
blue crosses) show good agreement. The scatter around the proto- 
cluster's redshift of approximately <x- = 0.3 is similar to the width 
chosen for the redshift cut in the LBG and BBG selection. Because 
of this good agreement we expect that the photometric redshift es- 
timates are adequate representations of the true redshifts. 

There are six additional objects with z spcc > 3, but only three 
of these are classified as LBGs. It thus seems that approximately 
half of the z ~ 3 galaxies have been missed have by the colour se- 
lection technique. We discuss this further in Sect. 14,51 Furthermore, 
the photometric redshift of the RG differs strongly from its spectro- 
scopic redshift. This is due to its proximity to a foreground galaxy 
which affects the photometry. 



4.5 Completeness: photometric redshift selection 

In addition to the selected protocluster galaxy candidates, we se- 
lected unclassified objects that have photometric redshifts in the 
range 2.7 < z p i lot < 3.5. The size of this sample gives an estimate 
of the number of z ~ 3 galaxies that are missed by the above- 
mentioned colour-selection techniques. 

Galaxies with 2.7 < z phot < 3.5 were selected from the R band 
and K s band detected catalogues. These objects must be located in 
the field-of-view covered by the 14 band, since the LBGs are se- 
lected from this band. They must also have NIR data in order to 
have the best constraints on the photometric redshift. All known 
LBGs, BBGs, ISAAC ghost images, spurious detections and dou- 
ble entries were removed and the appropriate magnitude cuts were 
applied. A sample of 96 objects was obtained of which 93 were 
detected in R band and 3 in either of the K s bands. 

Figure [8] shows the photometric redshift distribution of the R 
band detected galaxies. The distribution differs significantly from 
the Zphot distribution of LBGs as there is a large number of objects 
located at the lower edge of the redshift range. 

The majority of galaxies with 3.0 < z p h t < 3.3 is undetected 
in the 14 band and was assigned lower limits for their U± - V v 
colours. These 14 - V v lower limits are generally too blue for these 
objects to be classified as LBGs, but deeper U band data might lead 
to redder 14 - V v colours, shifting them into the LBG region. The 
inset in Fig.[8]shows the z p h 0l distribution of these possible LBGs. A 
peak at the redshift of the protocluster is seen. This peak accounts 
for a large fraction of the unclassified objects with 3.0 < z phot < 
3.3. There are 55 potential LBGs (pLBGs) that have lower limits 
on t4 _ V v that are too blue to be classified as LBGs. The LBG 
selection therefore misses up to 50 per cent of z ~ 3 galaxies. 

In the K s band only three additional objects are found to be 
missing from the BBG sample. The same redshift range yields 15 
unique BBGs. Based on this we estimate that the BBG selection 
yields ~ 80 per cent of all z ~ 3 K s detected galaxies. 

The positions of the objects in the final samples as well as the 
fields covered by the various instruments used are shown in Fig. [9] 
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Figure 9. Sky coordinates of the various populations. Large and small diamonds indicate the LBG and pLBGs, respectively. Blue crosses denote LAEs, green 
triangles [Ora] emitters and the squares denote the BBGs. A distinction is made between ISAAC and MOIRCS detected BBGs, with the former indicated by 
open red squares and the latter by filled yellow squares. The dotted lines indicate the outlines of the fields covered by the various instruments. The FORS data 
(Vfl) covers the entire figure and also includes all VIMOS bands (U V BV V R). Subsequently, the IRAC data covers the majority of the field with the exception 
of some of the corners. Going from larger to smaller fields we then have the field, the field covered by ACS, the two MOIRCS fields and finally, slightly 
below centre the area covered by ISAAC. 



5 SED FITTING 

Do the properties of the protocluster galaxies differ from those of 
field galaxies at the same redshift? To compare the protocluster 
galaxies to field galaxies we determine the galaxy properties by 
fitting the individual galaxy SEDs with population synthesis mod- 
els, and then compare the properties to equivalent samples detected 
in the field. 

The fitti ng of the galaxy SE Ds was done using the FAST code 
described in Kri ek et alj (2009). Due to its versatility and speed 
the use of the FAST code allows us to fit a large range of models 
and model parameters. Both the BC03 and the updated CB07 mod- 
els (Chariot & Bruzual, private communication) were used with a 
Salpeter IMF and Z =Z G . The inclusion of an improved treatment 
of the TP-AGB phase in the CB07 models results in a more reliable 
stellar mass estimate for populations older than 100 Myr. 

The free parameters in the fitting routine are the age, mass, 
SFH and the extinction by dust. We considered exponentially de- 
clining SFHs with decay times, r, ranging from 10 Myr to 10 Gyr 
with steps of 0.1 dex. The inclusion of SFHs with values of r much 
larger than the age of the Universe for our adopted cosmology al- 
lowed us to mimic a constant SFH. The age grid that we considered 
ranges from log(age/yr)=7 to 9.3 with age steps of 0.1 dex, where 
log(age/yr)=9.3 equals the age of the Universe at z ~ 3.13. The 
redshifts were fixed to the spectroscopic or photometric redshifts, 
depending on which was available. As the EAZY code is special- 
ized in determining redshifts and is more efficient in use, we have 
chosen to use the EAZY redshifts rather than including it as a free 



parameter during the SED fitting. This choice has little effect on the 
main results for the majority of the objects. 

The effect of internal d ust extinction the m odel SEDs was 
taken into account using the ICalzetti et al.1 (2000J) extinction law 
for values of A\ ranging from to 3 with steps of 0.1. The attenu- 
ation blueward of the Lyq line d ue to the IGM was included using 
the prescription of Madai] dl995h - 



6 RESULTS 

6.1 Number densities of galaxy populations 

6.1.1 Lyman Break Galaxy candidates 

To determine whether there is an overdensity of LBGs in the 0316 
field due to the presence of the protocluster, the number density of 
LBGs in the 03 16 field is compared to the ~ 0.3n° MUS YC ECDF- 
S blank field (see also Sect, \4.2j . Photometry for the control field 
was obtained as described in Sect. [3] 

The GOODS-S field, which is part of the ECDF-S field, ha s 
been shown to be underdense in DRGs ( I van Dokkum et all 2006), 
Comparison of the number density of DRGs and LBGs in 8 ECDF- 
S subfields reveals that the variation in LBG number density is 
lower than for DRGs. Furthermore, the central subfield, which has 
the lowest DRG number density, has a relatively high number den- 
sity of LBGs. We conclude that the ECDF-S is unlikely to be as 
underdense in LBGs as the GOODS-S field is in DRGs. 
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Figure 10. completeness-corrected cumulative number density of objects 
that satisfy the colour selection criterion as function of total R band magni- 
tude. Diamonds denote values for the 0316 field, whereas the squares indi- 
cate the results found for the MUSYC ECDF-S control field. Uncertainties 
are based on Poisson statistics. An excess of LBGs in the 0316 field is ap- 
parent. The inset shows the completeness of the fields as a function of mag- 
nitude with the full and dashed line denoting the 0316 and ECDF-S data, 
respectively. 

For the MUSYC ECDF-S data an LBG colour criterion was 
devised that is equivalent to the criterion used for the 0316 field: 

U-V> 1.85, 

V - R < 0.62, (3) 
U- V > 4.0 x(V-R) + 2.51, 
R < 26. 

Applying this to the MUSYC data yields a total number of 694 
LBG candidates. 

Figure[l0]shows the completeness-corrected cumulative num- 
ber density of LBGs for both the 0316 and the ECDF-S control 
field. Note that no information on the redshifts of the objects is 
used to select these LBG samples. There is an excess in the 0316 
field across the entire magnitude range. The 0316 field is a factor of 
1 .6 ± 0.3 denser than the ECDF-S field for galaxies with R < 25.5, 
with the l<x uncertainty based on Poisson statistics. Defining the 
galaxy surface overdensity as <5 g = «0316/«ecdf-s - 1, we find 
S g = 0.6 ± 0.3 in the 0316 field. This overdensity is not the result 
of a difference in general number counts, as the area-normalized 
completeness-corrected number counts of the two fields differ by 
less than 1.5 per cent. 

The significance of the LBG overdensity in 0316 is quantified 
by considering field-to-field variations of LBGs on the scale of the 
0316 field. The number density of LBGs in 125 subfields of the 
ECDF-S (each the same size as the 0316 field) is measured, and the 
resulting distribution displayed in Fig. QT| the arrow indicates the 
number density of LBGs in the 03 16 field. Only 1 in 125 fields has a 
number density that exceeds the 0316 density, thus the significance 
of the LBG overdensity in the 0316 field is at the 3<x level with 
respect to field-to-field variations. 

The surface overdensity of LBGs in the 0316 field is a lower 
limit to the volume overdensity of the protocluster in this field. The 
LBG selection criterion selects galaxies in a relatively large redshift 
range, so the LBG overdensity within the protocluster is larger than 
0.6 (see Sect.|7J. 

The outer VIMOS quadrants are used to assess the LBG num- 



Figure 11. Distribution of LBG number densities with R < 25.5 in the 
ECDF-S for 125 randomly chosen Keck-sized (40. In') fields. The arrow 
indicates the value found for the 0316 field. One in 125 ECDF-S subfields 
is found to host a larger number density than found in the 0316 field. 

ber density at greater distances from the RG. There is no i/ k cover- 
age for these outer fields, so the criterion given in Eq.|2]cannot be 
used. Instead, a similar criterion that uses the t/ v filter is applied. 
The number density in the outlying fields is approximately 10 per 
cent smaller than the central field. Therefore these outlying fields 
are overdense with respect to the ECDF-S, but to a lesser degree 
than the central 0316 field. This indicates that the protocluster ex- 
tends beyond the central field, out to at least 15' (~ 7 Mpc). This 
is larger than most esti mates of protocluster sizes which typicall y 
find sizes of 2-5 Mpc dlntema et alj|2006l IVenemans et al 1 (2007h . 
but often size estimates are limited by the field size. Nevertheless, 
the larger size may be a consequence of the possible superstructure 
in the field, as hinted at by the redshift distribution of the [O in] 
emitters. 

6.1.2 Balmer Break Galaxy candidates 

The study of K07 found that the 0316 field is overdense with 
respect to the GOODS-S field, with a pproximately 1.5—2 times 
more BBGs in the 0316 field. However. Ivan Dokkum et ai]( l2006h 
has pointed out that the GOODS-S field is underdense, contain- 
ing only 60 per cent of the number of objects found in th e larger 
MUSYC survey dGawiser et alll2006bl ; lOuadri et alj|2007t) . To as- 
certain whether there is an overdensity of BBGs in the 0316 field, 
we compare to the publically available MUSYC data described in 
IOuadrietal]d2007h . 

The four MUSYC fields cover a total area of ~400n'. Using 
these fields as control fields will yield better statistics compared to 
the GOODS-S, as well as decrease the influence of cosmic vari- 
ance. The control fields are denoted as HDFS1, HDFS2, 1030 and 
1255, respectively. Catalogues were constructed from the four im- 
ages following the procedure used for the 0316 field. Applying the 
BBG criterion and a K s < 23 cut yields 123, 116, 197 and 118 
BBGs for the HDFS1, HDFS2, 1030 and 1255 fields, respectively, 
or 554 BBGs in total. 

Fig. [T2] compares the completeness-corrected cumulative 
number density of BBGs in the 0316 field to the MUSYC control 
field. The small field-of-view ISAAC data suggests that the 0316 
field has an excess of bright BBGs (K s < 21.5) close to the radio 
galaxy, but the sample is small so the number statistics are poor. 
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Figure 12. Completeness-corrected cumulative number density as function 
of K s band magnitude for both populations of BBGs. ISAAC and MOIRCS 
detected BBGs are denoted by diamonds and asterisks, respectively. Uncer- 
tainties are obtained using Poisson statistics. The squares denote the mean 
cumulative number density of the four individual MUSYC fields. Neither of 
the 0316 BBG populations shows evidence for an overdensity in the 0316 
field. 
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Figure 13. Stellar mass distributions for each protocluster population. Ob- 
jects that are strongly affected by confusion in the IRAC images, or galaxies 
with M < 10 7 M© are not included. The total number of objects in each pop- 
ulation is shown in the upper left corner of each panel with the same colour 
coding as the respective histograms. The arrow indicates the best-fitting 
mass that is found for the LAE stack (Sect. l6~23V The hatched region indi- 
cates where mass estimates become highly uncertain due to faintness in the 
restframe optical. 



At fainter magnitudes there is no overdensity. Since the MOIRCS 
field shows no indication of an excess at bright magnitudes, we 
conclude that there is no evidence for an overdensity of BBGs in 
the 0316 field. This is in contradiction with the conclusion of K07. 

The lack of an overdensity seems to contradict the presence of 
a peak in the photometric redshift distribution at the protocluster's 
redshift (Fig.[%J». Nine ISAAC selected BBGs lie in the range 2.7 < 
Zphot < 3.5. Only four of these have K\ < 23, so more than half of 
the BBGs at the protocluster's redshift are too faint for inclusion 
in Fig. [12] It is possible most of the protocluster BBGs are faint, 
and would only be revealed as an overdensity at faint magnitudes, 
however the control fields are not deep enough to warrant a proper 
comparison. 

6.2 Properties of protocluster galaxy candidates 

Below we present the results of the SED fitting for all protocluster 
populations. The LAEs are treated separately in Sect. 16.2.51 be- 
cause they are generally very faint across the entire wavelength 
range. This results in poorly constrained SED fits and therefore a 
different approach was taken for this population. 

6.2.1 Mass 

The mass distributions of all protocluster candidate populations are 
shown in Fig.[T3] The estimate of the galaxy's stellar mass is domi- 
nated by the flux in the IRAC bands, so all objects that are strongly 
contaminated by neighbouring sources have been excluded from 
the samples. Table[2]lists the median mass values for each popula- 
tion. 

The upper right panel of Fig. [T3] shows that most LBGs have 
masses between 10 9 and 10" M Q and the median stellar mass is 
a few times 10 9 M Q . This mass increases to ~ 2xl0 10 M Q when 
limiting the sample to LBGs with R < 25.5. The median mass of 
the pLBG population is approximately a factor 2 lower than the 
LBGs, and a KS test shows that the mass distributions differ at the 
2tr level. Thus the z ~ 3 galaxies that are not included in the LBG 



sample due to insufficient depth of the U band are generally less 
massive than the LBG galaxies. 

A large fraction of the [Om] emitters is detected in the [3.6] 
and [4.5] bands, whilst a few are also detected in the [5.8] and [8.0] 
bands. Six objects (the RG and source IDs 2, 6, 10, 11 and 12 in 
M08) are strongly contaminated by nearby bright objects. The re- 
maining 7 objects are all detected in the [3.6] and [4.5] bands. This 
implies that the [Om] emitters have significant stellar masses. The 
median stellar mass of the [Om] emitters is larger than the median 
LBG mass, however, this median mass is based on a sample of only 
7 objects. 

[Om] emitters are selected as high equivalent width objects in 
the observed K s band, thus they are selected to have low K s band 
magnitudes relative to the [Om] narrowbanc^f]. Hence it is surpris- 
ing that such a selection criterion identifies massive galaxies. 

Not all [Om] emitters have been spectroscopically confirmed, 
so it is possible the sample contains low redshift interlopers. The 
prime suspects for interlopers would be Ha emitters at z ~ 2.15. 
Three out of the 7 [Om] emitters with reliable mass estimates 
are not spectroscopically confirmed and two of these have z p i lot ~ 
2.17. Moreover, all spectroscopically confirmed [Om] emitters have 
Zphot > 3. Thus it is possible that these two objects are interlopers. 
Removing these objects from the sample leaves three galaxies with 
M =6-8xl0 9 M G and two galaxies with ~ 10" M G . 

The BBGs have masses of a few times 10 10 M Q with some 
exceeding 10" M Q . The BBGs comprise the population with the 
highest masses of the protocluster candidate populations. Since the 
BBGs are K s band selected it is expected that they will have high 
stellar masses. The ISAAC-selected BBGs are on average less mas- 
sive than the MOIRCS-selected BBGs, but this is only due to the 
larger depth of the ISAAC data. When the same magnitude cut 
is applied to both datasets the mass of the ISAAC-selected BBGs 
agrees with that of the MOIRCS-selected BBGs. 



15 At low redshifts the [Om] line is often found to be an indicator of AGN 
activity. Therefore part of the IRAC flux may be caused by heated dust close 
to an AGN. 
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6.2.2 Age and extinction 

The SED fitting procedure results in degeneracies between the best- 
fit age, extinction and SFH. Because of the large uncertainties in 
the best-fit ages of the LBGs, we do not discuss them further in this 
work. The ages of the BBGs are easier to constrain, because they 
were selected to have strong Balmer breaks. The median age at 
which the BBGs started forming stars is — 1 Gyr ago (see Table[2}. 



6.2.3 UV slope 

Assuming the restframe UV continuum has the form f A = CAP, we 
calculated the UV slope B of the candidate protocluster members. 
Either the R and / or r 62 s and / 814 bands were used, depending on 
whether ACS coverage is available. At z ~ 3 these bands corre- 
spond to restframe 1500 A and 2000 A, respectively. The UV slope 
was calculated as 

logjfo -l°g/i,i 



l0g/le: 



log^e 



(4) 



with /lcfr,R/i the effective wavelength of the R or / filters, respec- 
tively. The uncertainty on B was obtained by varying the flux by an 
amount drawn randomly from the normal photometric error distri- 
bution. 

The median B for the LBG population B is -1.52, but the bright 
subset (R < 25.5) has a redder UV slope of -1.38. This reddening 
of 8 w ith increasing brightness was also found bv lBouwens et al.l 
(2009). An unobscured, young population has a UV slope ranging 
from -2 to -2.5. Thus the LBGs are slightly attenuated by dust. 
The pLBGs are on average bluer than the LBGs, indicating there is 
less dust in these objects. 

The [Om] emitters have similar UV colours to the LBGs. 
Whereas the BBGs have significantly redder colours than the LBGs 
so they are either more heavily obscured or older than the LBGs. 



6.2.4 Star formation rate 

The star formation rates (SFRs) of the candidate protocluster galax- 
ies were esti mated from the res t-frame 1500 A luminosity using the 
relation from lKermicuttl i 19981) 



SFR vw [M Q yr- 



(5) 



8 x 10 27 ' 

This relation is valid for a Salpeter IMF and does not take the pres- 
ence of dust into account. For simplicity we assume that all objects 
are located at the redshift of the protocluster. Furthermore, this rela- 
tion is only valid for starbursts that have reached an equilibrium be- 
tween the number of stars forming and evolving off of the main se- 
quence. The SFRs of extremely young starburst galaxies that have 
not yet reached this equilibrium will be underestimated. 

The SFRs of the candidate protocluster galaxies are shown 
in Fig. [14] Most galaxies have SFRs less than 10 M Q yr" 1 , ir- 
respective of galaxy population, with a small fraction of outliers 
having SFRs up to 40 M Q yr" 1 . Both the median dust-uncorrected 
and dust-corrected SFRs for each population are listed in Table [2] 
The extinction was derived from th e UV slopes using E(B - V) = 
03 - yg )/8.067 jMeurer et alj|l995h and A(1500) = 4.39£(S - V) 
( Calzetti et afl l2000l) . Here B is the UV slope for an unattenuated 
ionizing population of stars, assumed to be -2.5. This method of 
estimating SFRs gives good agreemen t with SFRs derived from ra- 
dio observations for galaxies at z ~ 2 jPannella et all2009I) . 

The LBGs have a relatively small amount o f dust extinc- 
tion with a median A(1500) = 0.53 (Ay = 0.21). Ishaplev et al.l 
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Figure 14. Extinction uncorrected star formation rates for the various popu- 
lations obtained from the restframe 1500 A luminosity. The arrow indicates 
the SFR obtained for the LAE stack. For simplicity, we have assumed that 
all LBGs and BBGs with 2.7 < z < 3.5 are located at the protocluster's 
redshift. 



d200lh found generally larger dust extinction with a median value 
of A v ~ 0.6. We revisit this difference in a more thorough compari- 
son in Sect. |7.2.2l No significant difference was found between the 
uncorrected SFRs of the LBGs and pLBGs, but the dust-corrected 
SFRs of the pLBGs are systematically lower, consistent with the 
UV slopes of the pLBGs being systematically bluer. 

The slope of the UV continuum can become redder due to 
both aging of the stellar populations, and attenuation by dust. Thus 
it is ambiguous to use the UV slope as an indicator of extinction 
for BBGs. Therefore the amount of extinction was measured from 
both the UV slopes and the value for A{ 1 500) derived from the SED 
fitting. Both values of the dust-corrected SFR are listed in Table[2] 

The estimated SFRs of the MOIRCS BBGs are systematically 
larger than those of the ISAAC BBGs. However, the MOIRCS data 
are shallower than the ISAAC data and therefore include objects 
with systematically larger mass. The observed difference in SFRs 
between MOIRCS and ISAAC objects is co nsistent with the in- 
crease of SF R with increasing mass at z < 1 dNoeske et ai]|2007l) 
and at z ~ 2 jDaddi et alj|2007l) . 



6.2.5 Lya emitters 

In the sample of LAEs, 27 are not detected in any IRAC band, 
indicating that the LAEs have small stellar masses. Figure[T3lshows 
that most LAEs have M < 10 9 M Q . The properties of these IRAC- 
undetected LAEs were investigated by mean stacking 16 LAEs that 
lie within the field-of-view of the U V BVRI and MOIRCS datasets 
and are not contaminated by bright neighbours in the IRAC data. 

The photometry and the best-fitting SED for the stacked data 
are shown in Fig. [15] The stack reveals robust detections in the 
BVRI bands, but not in bands redward of /. The derived stellar 
mass is approximately 10 8 M Q and this is marked in Fig. 1131 with 
an arrow. It is much smaller than the average LBG stellar mass. 
The best-fit age is 10 Myr, which is the lowest age allowed by 
the SED fitting process. The amount of dust extinction is uncon- 
strained from the fit to the SED, but the UV slope determined from 
the rs25 and /si4 bands is = -2.4. Thus the average protocluster 
LAE contains little dust. These results are in agreement with pre- 
vious studies that show that LAEs are in general young objects ex- 
periencing their first burst of star formation (e.g.lHu & McMahonl 
1996; Mal hotra & Rhoadsll2002l ; lTapken et al.ll2004h . 
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Table 2. Upper part: median values of various properties for the samples discussed in this work. Uncertainties here are mean absolute deviations on the median 
divided by the square root of the number of objects. Lower part: Properties for the stacked and bright, individual LAEs. 



Object log(M)cB07 log(M) Ba)3 SFR (M Q yr" 1 ) SFR C01T (M G yr" 1 ) log(Age) CB 07 P 



LBG 9.92 10.1 7.0 11.3 - a -1.52±0.07 

LBGfi <2 5.5 10.16 10.29 11.4 23.1 - a -1.38±0.09 

pLBG 9.65 9.79 5.7 9.1 - a -1.72±0.08 

pLBG« <25 .5 9.79 9.95 8.7 13.9 - a -1.54+0.15 

Ora b 10.51 10.62 6.3 9.5 - a -1.50±0.13 

BBGi 10.55 10.77 2.2 C 5.1/6.9 cd 9 -1.12±0.21 c 

BBG M 11.0 11.15 6.6 C 15.0/6.8 cd 8.9 -0.91±0.19 c 



LAE stack 8.1+°" 0.9 1.0 7.0 +1 ° -2.40±0.43 

LAE, 5 ,8 9.97+gJi, lO.O?^]^ 15.7 20.4 8.8+!M -1.97±0.32 

LAEi 867 1 1 -08+^ °| 11.22™ 40.2 64.2 8.8+° : ' -1.56±0.73 

LAE 310 1 10 - 48 !o68 10 - 74 !o94 n - 3 27 - 6 8 - 8 -l7 -0.7±1.8 



a Individual ages are unconstrained 

b Values for the SFRs and /? taken from M08 

c Only objects with 2.7 < z < 3.5 are taken into consideration 

d First value calculated using the UV slope and second value calculated using the SED fit results 
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Figure 15. Photometry of a stack of 16 LAEs and the corresponding best- 
fitting SED. The diamonds indicate the photometry as computed from the 
best fit SED. Best-fitting parameters are listed in the bottom right corner, 
lcr upper limits are used for [5.8] and [8.0] as the stack yields negative flux 
in these bands. The best-fitting model is very young and has a low mass, 
but a significant amount of dust extinction. Due to the large uncertainties on 
the JHK and IRAC photometry a model SED without dust would lead to an 
equally acceptable fit. 



Only 4 of 32 LAEs have robust detections in [3.6] and [4.5], 
whilst 1 LAE has a marginal detection. The four objects with sig- 
nificant IRAC detections include the RG. However, since the RG 
is located close on the sky to a foreground galaxy, the photometry 
in all bands is compro mised. It is therefor e not included in further 
analysis (but see also lSevmour et al.l2007l who estimate its mass to 
be 1.6x10" M Q based on a Kroupa IMF). The remaining 4 LAEs 
are discussed below using the IDs designated in V05. 

1518: The detection of this object in the [3.6] and [4.5] bands 
is marginal. However, it has significant detections in both / and K[ 
(6cr and 7.5cr, respectively). The photometry and the best-fit SED 
are shown in the upper left panel of Fig. [T6] The galaxy is best 
fit by a model with continuous SFR and an age of 600 Myr. No 
dust attenuation is required. The best-fit stellar mass is ~ 10'° M Q , 



making it almost two orders of magnitude more massive than the 
average LAE. 

1867: Excluding the RG, this object is the second brightest 
LAE in the sample. It has a clear IRAC detection in all four bands 
and is also bright in the restframe UV/optical continuum. The pho- 
tometric redshift of this object is 3.15, which is consistent with the 
spectroscopic redshift of 3.134. The galaxy is of similar age to ob- 
ject #1518 with a best fit age of 600 Myr and r ~ 200 Myr. It 
contains no dust and it is very massive, having a stellar mass of ap- 
proximately 10" M Q . This galaxy is also selected as an LBG and a 
BBG It has a clumpy morphology as seen in ACS data (V05, Ven- 
emans et al. in prep.) and is located 50" (or 350 kpc in projection) 
from the RG. 

2487: This object hosts an AGN, as evidenced by the broad- 
ness o f the Lyq line and th e presence of a broad CIV emission line 
(V05. lLeFevre et al.ll9 96). Since the SED fitting does not account 
for an AGN contribution, the best fit parameters are not valid and 
are not discussed further. 

3101: This LAE has similar properties to #1867, but it is less 
massive by a factor 3. Because of the absence of ISAAC data, the 
best-fit properties are poorly constrained. 

An overview of the properties including lcr uncertainty inter- 
vals can be found in Table [2] Note that the two objects for which 
only MOIRCS NIR data is available (#2487 and #3101) have sig- 
nificantly larger uncertainties on the best-fitting age. This reflects 
the aforementioned need for deep NIR observations in order to con- 
strain the age. 



6.3 Dependance of galaxy properties on location within the 
protocluster 

An important question in studying the history of cluster evolution 
is whether there are systematic variations in the properties of the 
protocluster members that depend on their location within the pro- 
tocluster, e.g. distance from the radio galaxy. 

In Fig.[l7]the stellar mass of the LBGs is plotted as a function 
of distance from the RG. The main panel of Fig.[T7]shows the area- 
normalized integrated mass M illt within annuli of 50", whilst the 
inset shows the individual masses and the mean mass per annulus. 
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Figure 16. Photometry and SED fits for the four LAEs that were detected in the [3.6] and [4.5] bands. Best-fitting parameters are shown in the bottom right 
corner of each panel. Diamonds denote magnitudes as calculated from the best-fitting SED. ID numbers correspond to the ones used in V05. The first two 
objects fall inside the limited ISAAC field-of-view and are shown with the ISAAC / and K s data for the restframe optical. For the remaining two objects no 
ISAAC data was available and the shallower MOIRCS JHK S data is used. Three out of four SED fits return ages of 500 Myr and higher and all masses are 
approximately larger than 5 X 10 9 . Non-detections are replaced with lcr upper limits. 



Two massive LBGs near the RG result in a high M ml in the inner- 
most annulus. M mx then sharply decreases out to 150", but slightly 
increases at greater distances. 

The inset in Fig.ll7lshows a large scatter in the masses of the 
LBGs at all radii, which is expected with such a large number of 
interlopers. This large scatter in combination with the small sample 
size could also cause the upturn at distances greater than 150". 

The two most massive LBGs are located near the RG: one 
is a confirmed LAE protocluster member, whereas the other is an 
[Om] emitter with z p hoi = 3.23. Both galaxies have stellar masses of 
approximately 10" M G making them almost as massive as the RG 
itself (2.5-3x10" M Q for a Salpeter IMF). This is different from 
what is found for radio-selected protoclusters at z =4-5, where the 
RG is approximately an order of mag nitude more mass ive than the 
brightest of the surrounding LBGs dOverzier et al.|[200gh . 

Figure[T8"lshows the area-normalized integrated SFR in annuli 
of 50". A similar trend is observed for the SFR as for M mt : galax- 
ies close to the RG have larger SFRs, although again the scatter is 
large. At small distances (< 100") all objects are forming stars at 
rates of >15 M Q yr -1 , whereas at larger distances the majority of 
the objects are forming stars at ~5 M Q yr~' . That galaxies with the 
highest SFR are located in the centre of the protocluster is in stark 



contrast with the situation in galaxy cl usters at z < 1 where the 
SFR d rops towards the cluster core fe.g.lBalogh et alJll998l,ll999l : 



lHashimoto et al]l998l ; |Postman et alfeOOlklLewis et al.ll2002l) . 

The above figures show that the RG is surrounded by some 
of the most massive and active galaxies in the protocluster. How 
unique is the position of the RG in the protocluster? The above pro- 
cedure was repeated for 200 random locations in the 0316 field, and 
the resulting mass and SFR distributions for the innermost annulus 
are shown in Fig. [19] There are few locations that yield similar or 
larger M mt or total area-normalized SFR, and these are all located 
close to the RG. Thus the immediate surroundings of the RG are 
unique within the 0316 field. The radio galaxy likely lies in the 
core of the forming protocluster and the protocluster core galaxies 
are affected by their environment. 

Note that there are several effects that will weaken such ob- 
served spatial trends. First, since the photometric redshifts alone 
cannot establish protocluster membership, there will be contamina- 
tion by interlopers in such an analysis. The observed LBG surface 
overdensity of 0.6 implies that slightly less than 2 out of every 3 
LBGs is an interloper. Sec ondly, many protocluster systems show 
filamentary structure s (e.g. lOuchi et al . 2005; Matsud a~et alj|2005l : 
lOverzier et al J 20081) . hence the assumption of a circular symmetry 
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Figure 19. Left panel: Distribution of total area-normalized mass within a 50" radius for 200 random locations across the 0316 field. Right panel: Same as 
left panel but for SFR. Arrows indicate the total area-normalized mass and SFR in a 50" annulus when the RG is taken as the centre of the protocluster. 
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Figure 17. Stellar mass of all LBGs as function of distance to the RG. 
Masses of individual objects within annuli of 50" width are summed and 
normalized to the area of the annuli. The inset shows the individual values 
as function of distance to the RG and mean values for each of the annuli. 
The individual values show a weak trend, but the integrated mass shows a 
sharp decrease as the distance to the RG increases. 



Figure 18. Total area-normalized SFRs for annuli of 50" width. SFRs are 
dust-uncorrected. The inset shows the individual values and the mean val- 
ues for bins of 50". Vertical error bars are mean absolute deviations. The 
seemingly lower limit to the SFR present in the inset is due to the magni- 
tude cut at R = 26 used to select LBGs. Objects close to the RG are forming 
stars more rapidly. 



used above will also weaken the signal. Thirdly, the RG itself is not 
included in this analysis, and its large mass (> 10" M Q ) and SFR 
(110 M Q yr~' without dust correction) will lead to an even stronger 
correlation if included. Thus the radial dependence is likely to be 
stronger than observed. 

We conclude from the observed radial trends in stellar mass 
and SFR that galaxy evolution is strongly influenced by the proto- 
cluster environment, with larger SFRs and stellar masses occurring 
in the protocluster core. 



7 DISCUSSION 

7.1 Surface and volume overdensities 

7.1.1 Lyman Break galaxies 

The measured LBG surface overdensity of the 0316 field is 0.6 ± 
0.3. Because the uncertainties of the photometric redshifts (<r z /(l + 
z)) range between 0.02 and 0.35, the derived surface density will be 



diluted by the fore- and background objects that are in the sample 
as well. The surface overdensity of 0.6+0.3 that is found is thus 
a lower limit to the volume overdensity of the LBGs. How strong 
must the volume overdensity of LBGs be to produce a surface over- 
density of 0.6? 

V05 showed that the FWHM of LAE redshifts in the 0316 pro- 
tocluster is Az ~ 0.025. LAEs are in general significantly less mas- 
sive than LBGs so they may relax within the protocluster potential 
more easily. We assumed that the protocluster spans a line-of-sigh t 
distance of Az = 0.04 (as found in the study of lSteidel et al.l (1998) 
for LBGs in the SSA22 protocluster at z = 3.09). In the field-of- 
view covered by the band, the protocluster spans a volume of 
~ 5.6 x 10 3 Mpc 3 . 

The LBG selection criterion selects galaxies in the redshift 
interval 2.9 < z < 3.4, as shown in the inset in Fig. [3] Although 
the spread in the z p hot distribution exceeds this interval, this can 
be attributed to the uncertainties in the z p hot estimates. There are 
448 LBGs with R < 25.5 in the ECDF-S and these span a volume 
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of 1.8xl0 6 Mpc 3 . This implies that there should be approximately 
16-17 'field' LBGs in the 0316 field. 

To obtain a surface overdensity of <5 g = 0.6 an additional 10- 
1 1 LBGs are needed in the protocluster volume. The protocluster 
volume can host approximately 1.3 'field' LBGs. Hence the proto- 
cluster is a factor ~ 8 ± 4 denser than the field and the volume over- 
density S g y is 7 ± 4, with the l<r uncertainties calculated by prop- 
agation of the uncertainty on the surface overdensity. This is larger 
than the overdensities observed in the LAEs and [Om] emitters, but 
this volume overdensity is little more than an order of magnitude 
estimate due to the many uncertainties in the assumptions that were 
made. 

Furthermore, based on the presence of three [Om] emitters at 
z = 3.1, M08 speculated the 0316 protocluster is part of a larger 
60 comoving Mpc superstructure. This superstructure would ex- 
tend beyond the redshift range probed by the hya NB filter. If the 
protocluster is part of a larger superstructure, the protocluster vol- 
ume assumed above would result in an overestimate of the LBG 
overdensity. 

The mass of the protocluster was computed using the comov- 
ing volume V occupied by the protocluster, the mean density of 
the Universe p and the matter overdensity 6 m which relates to the 
galaxy volume overdensity S g y as S m = 6 s y/b with b being the bias 
parameter. Using V - 5.6 x 10 3 Mpc, p = 1.9 x 10'° M a Mpc' 3 
and b =1-5 (e.g , . lAdelbereer et alj[l998l : iGiavalisco et"atlll998l ; 
lOuchi et al.ll2004! ; lLee et alj|2006h the protocluster mass is M =2- 
12xl0 14 M Q . The range was determined by the combination of the 
lcr uncertainty and the range in the bias value. This protocluster 
mass is consistent with V05 and is similar to the mass of the Virgo 
cluster. Since it was shown in Sect. 16.1.11 that the overdensity is 
spread over multiple VIMOS quadrants, the protocluster volume is 
likely underestimated. Therefore, the cluster mass should be con- 
sidered a lower limit. 



7. 1.2 Balmer Break galaxies 

We found no surface overdensity of BBGs in the 0316 field based 
on the number densities. But would we expect to detect a sur- 
face overdensity in the BBGs assuming the protocluster has the 
same volume overdensity of BBGs as it has of LBGs? Or will the 
overdensity be smoothed out due to the large redshift range that is 
probed by the BBG selection criterion? 

The 554 BBGs detected in the four MUSYC fields are as- 
sumed to lie between 2 < z < 3.5 and are homogeneously spread 
across this redshift interval. Thus the volume spanned by this popu- 
lation is 2.2x 10 6 Mpc 3 . This implies that there should be 7-8 'field' 
BBGs with A: s < 23 in the ISAAC field-of-view. The protoclus- 
ter occupies 755 Mpc 3 (within the ISAAC field-of-view), and this 
volume contains on average 0.2 'field' BBG. So if the BBGs are 
overdense by the same factor as the LBGs (volume overdensity of 
7), the protocluster will host ~ 1.4 BBG. This results in a surface 
number density in the ISAAC field that is 1.2 times as high as the 
field density and thus a surface overdensity S g = 0.2. This is consis- 
tent with the lcr uncertainties on the cumulative number densities 
at K s = 23 as shown in Fig. [12] We conclude that a volume over- 
density of 7 cannot be ruled out due to the small sample size. 

It is doubtful whether systematic searches for BBG overden- 
sities at z ~ 3 have any chance of succeeding considering the range 
in redshift that is probed. To evaluate this the amplitude of field- 
to-field variations for an ISAAC-sized field is assessed. A similar 
method as described in Sect. 16.1. His emplo yed. As reference the 
UDS catalogue from [Williams et al. (2009) is used. This field is 
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Figure 20. Distribution of BBG number densities with A^ s < 23 in the 
UDS field IWiUiams et aj|2009li for 1000 randomly chosen ISAAC-sized 
(5.4n' fields. The arrow indicates the respective number density found for 
the 0316 field. To obtain a 3cr overdensity signal with respect to the field- 
to-field variations a number density of at least 3.5/n'is needed. 



complete down to K s = 23 and covers a total of ~ 0.8n°. The BBG 
number density was evaluated in 1000 ISAAC-sized fields placed 
randomly across the UDS field. The resulting distribution is shown 
in Fig. [20] The arrow indicates the BBG density found in the 0316 
field. It is clear that this result is not significant with respect to the 
field-to-field variations. In fact, a number density of ~ 3.5 per □' is 
needed to obtain a surface overdensity that is significant at the 3cr 
level. This would imply a surface overdensity c5 e ~ 1.2, equiva- 
lent to a volume overdensity of ~ 45. Such a strong overdensity is 
highly unlikely to exist over a large field-of-view. We conclude that 
selection criteria that span a considerable redshift range are poorly 
suited for the detection of overdensity signals around HzRGs. 



7.2 Field and cluster ensemble properties 

Does the evolution of galaxies in crowded environments in the 
Early Universe, such as protoclusters, differ significantly from 
galaxy evolution in the field? At z = 2.3. ISteidel et alj ( 120051) found 
that both ages and masses of galaxies in a protocluster are ~ 2 times 
larger than similarly selected galaxies in the field. A similar study 
at z =4-5 has shown no significant differe nce between protoclus- 
ter and field galaxies dOverzier et al.l l2009). This study can give a 
better insight in the situation at z ~ 3. 

We compare galaxy masses, SFRs and UV slopes for the var- 
ious protocluster populations with those of equivalent populations 
in the field. There is no sample of field [Om] emitters in the liter- 
ature, hence this population will not be discussed. We also refrain 
from comparing the 0316 BBGs to field BBGs as there is no strong 
evidence for an overdensity in the 0316 field. The SED fit results 
given here have been obtained using BC03 models with a Salpeter 
IMF to facilitate comparison with previous studies. 



7.2.1 Lya emitters 

Studies of L AEs at z ~ 3 find stellar masses ranging be t ween 10 8 
and 1 9 Mo dGawiser et ai]|2006al ; iNilsson etal]|2007l ; lLai et al.1 
2008) which is fully consistent with the results presented in this 
study. 

We also presented SEDs for a small population of massive 
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LAEs that do not fit the picture of galaxies that experience their 
first burst of star formation. This subpopulat ion of massive L AEs is 
found frequently in recent studies. iFinkelstein et al, (2009) found 
that two objects in their sample of 14 LAEs at z = 4.4 were best fit 
by a model SED with an age of 800 Myr and a mass of 6.5x 10 9 M Q . 
This fraction of 2 out of 14 is similar to t he fraction of massive ob- 
jects found in this study. lLai et alj f2008h find a fraction of massive 
LAEs of approximately 30 per cent, but they quote a 2<x detec- 
tion limit of 25.2 for [3.6]. T his is approximately h alf a magni- 
tude deeper than the 0316 data. lPentericci et al] d2009) investigated 
Lycf-emitting LBGs, finding a small fraction of objects that are sig- 
nificantly older than the average LAE. In summary, there is no ev- 
idence that the LAEs in the protocluster are significantly different 
from the field population. 



7.2.2 Lyman break galaxies 

Taking all LBGs with/? < 25.5, the median mass of the 0316 LBGs 
is ~ 1.6 times more massive than that of 'field' LBG s, which have 
a median mass of 1.2 x 10 10 M Q Jshaplev et al]|200ll) . In addition, 
LBGs close to the RG are more massive than LBGs at larger radii, 
indicating protocluster environment promotes the increased stel- 
lar mass of its members. However, the BC03 populat ion synthesis 
model s were not available at the time of the study of Shapl ev et al.l 
d200l[). so a fully consist ent comparison can not be made. 

Lehmer et al J ( 120090 find that LBGs suspected to be in the pro- 



tocluster SSA22 at z = 3.09 have on average H band luminosities 
that are 1.2-1.8 times larger than LBGs that are unrelated to the 
protocluster. Since the H band can be interpreted as a proxy for 
stellar mass at z ~ 3 this also indicates that the masses of the pro- 
tocluster LBGs are larger by the same factor. This factor of 1.2-1.8 
is similar to the difference in median masses found in this study. 

The 0316 LBGs have dust-uncorrected SFRs of ~7- 
11 M G yr~' (depending on whether the whole sample is consid- 
er ed or only the objec ts that have R < 25.5). The R data presented 
in lSteideiet al. (2003) yields a median dust-uncorrected SFR for 
field LBGs of ~ 10 M Q yr -1 . This is consistent with the median 
dust-uncorrected SFR found for the 0316 LBGs. If the SFRs are 
corrected for dust extinction, the SFRs increase on average by a 
factor of 2. This small increase in SF R indicates the prese nce of 
small amounts of dust in the LBGs. Ishaplev et al.l ( feOOll) used 
SED fitting to find that field LBGs have a median extinction of 
E(B - V) = 0.155 (A v = 0.62). This is larger than found in the 
03 16 protocluster as half of the 03 16 LBGs are best fit by A v < 0. 1 . 
However, A v is poorly constrained, with 31 out of 48 LBGs having 
lcr uncertainties that are consistent with A v = 0.6. We can there- 
fore not conclude whether the 0316 LBGs are less dusty than field 
LBGs. 

iMeurer et alj d 19991) co mposed a sample of z ~ 3 LBGs 
from the HDF V2.0 catalog ^Williams et all 1 19961) using a UVI 
filter combination. This filter set is similar to the set used in this 
study. They find a median UV slope (ft) for field LBGs of approx- 
imately -1.6, similar to 0316 LBGs. We thus find no evidence for 
a diffe rence in the UV slo pes between field and cluster LBGs at 
z ~ 3. IShaplev et alj {2003) find significantly redder slopes for field 
LBGs, but the filter set used in their work differs significantly from 
this study, so it does not offer a fair comparison. 

We find that the 0316 LBGs have stellar masses, SFRs and 
UV slopes that are consistent with those of field LBGs. We thus 
conclude that at z ~ 3 the ensemble properties of all possible pro- 
tocluster galaxies is not significantly different from what is found 
in a low-density environment. 



7.3 The future of MRC 0316-257 

We have found that the 0316 field contains an overdensity of star 
forming galaxies. How will the protocluster LBGs evolve and at 
what epoch will a quiescent population emerge? 

Using a simple 'back of the envelope' calculation we estimate 
how the protocluster will appear at later times. Th e peak of the 
cosmic SFH lies below z ~ 3 (e.g. Steid el etalJl999T) and therefore 
new LBGs will emerge at later times. These 'new' objects are not 
taken into account and we only evolve the LBGs that are currently 
identified. The sample is therefore static and no new objects are 
introduced. 

The SED fitting yields both a best-fitting age and SFH, so it is 
possible to evolve the LBGs according to this information. We find, 
however, that both these quantities are poorly constrained. Hence, 
we make the following assumptions. Since the distribution of SFHs 
peaks at approximately r = 10 8 yr, we assume that the typical LBG 
has an SFH of 100-500 Myr. An age of 50-500 Myr is assumed, 
in accordance with earlier s tudies that state 100 Myr to be the typ- 
ical age for LBGs at z ~ 3 ( Sha plev et"al]|200ll) . Finally, the initial 
SFR is assumed to be 40-60 M Q yr~' . These SFRs are to be inter- 
preted as being intri nsic SFRs. They are based on the results from 
IShaplevet allfcoOll) and as shown in Sect. |7.2| are consistent with 
the results for the 0316 LBGs. For simplicity the influence of dust 
extinction will not be taken into account in this scenario. 

If a quiescent galaxy is defined as having an SFR < 2 M G yr~' 
we find that the majority of a population of LBGs will have reached 
quiescence after a period of approximately 1 Gyr. This corresponds 
to z ~ 2. Only the youngest LBGs having t = 500 Myr will last 
until z ~ 1.5 before they reach quiescence. Therefore, the 0316 
protocluster should harbour a substantial 'dead' population after 
1 Gyr of evolution. 

These predictions can be compared to observations of proto- 
cluster structures at lower redshifts. One of the best studied pro- 
toclusters is the overdense regi on surroundi ng the z = 2.16 radio 
galaxy PKS 1 138-262 (1 138) iMilev et al.l2006h . This radio galaxy 
has a intricate, clumpy morphology and seems to be in an advanced 
stage of merging with several smaller galaxies. It is surrounded by 
a diffuse UV halo which is e vidence for in-situ star formation out- 
side the main radio galaxy dHatch etalj|2008l) . Furthermore, the 
field has been shown to harbour overdensities in both Lyor and 
Ha emitting galaxies, extremely red objects (EROs), X-ray emit- 
ters an d sub-mm bright galax ies I Pentericci et all 200(31; Kurk et all 
2004bllal; Istevens et al. 2003]; ICroft et alj|2005l; Izirm et alj|2008l; 



Kodama et alj|2007h . iKurk et alj (2004b) has also shown that the 
EROs and Ha emitting galaxies are located closer to the radio 
galaxy than the bluer star forming galaxies. The 1138 protocluster 
is therefore a prime example of a z ~ galaxy cluster progenitor. 

IZirmet "all (2008) have shown that the 1138 system contains 
red galaxies, some of which might be quiescent. However, the red 
sequence is still in the stages of formation and has n ot been properly 
establi shed yet. This is in agreement with s t udies by Blakesle e et alj 
(2006) and lvan Dokkum & van der Mareil d2007l) that state that the 
formation redshift of the red sequence lies in the range z =2-2.5. 
This roughly agrees with our toy model for the 0316 protocluster. 
A significant number of red galaxies should be in place as early as 
z = 2.5. This is earlier than what most studies predict but we are 
limited by the simplicity of this calculation. If the SFH is episodic 
and burst-like the time it takes for an LBG to become quiescent 
increases, so we expect the red sequence will form at a later time 
than estimated by this toy model. 

We predict that the population of galaxies currently residing 
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in the 0316 protocluster will resemble the 1138 protocluster after 
1 Gyr of evolution and 0316, like the 1 138 protocluster, will evolve 
into a rich galaxy cluster at z < 1 . 



8 SUMMARY AND CONCLUSIONS 

We have presented a comprehensive study of several galaxy pop- 
ulations in the protocluster associated with MRC 0316-257 at z = 
3.13, using restframe FUV to optical images. In addition to study- 
ing the LAEs and [Om] emitters found in previous studies, we iden- 
tify samples of LBG and BBG candidates whose properties provide 
additional evidence for the presence of a protocluster. 

(i) The cumulative number density of the LBG candidates in the 
0316 field is a factor of 1 .6+0.3 larger than for comparable non- 
protocluster fields at similar redshifts, indicating a surface over- 
density of 0.6±0.3 for LBGs in the protocluster. The surface over- 
density is significant at the 3<x level with respect to the field-to-field 
variations. This measured LBG surface overdensity strengthens the 
conclusion of V05 that a protocluster surrounds MRC 0316-257. 
Using estimates of the protocluster size, the measured surface over- 
density gives an LBG volume density that is 8+4 larger than that 
of the field. Such an overdensity implies a minimum mass for the 
protocluster of 2-12xl0 14 M G . 

(ii) The redshift distribution of the BBG candidates shows a 
peak at the redshift of the radio galaxy. This is further evidence that 
there is an overdensity of galaxies and a protocluster around MRC 
0316-257. There is no significant surface overdensity of BBGs, but 
this is not surprising due to the small sample size. A volume over- 
density of BBGs comparable with the measured volume density of 
LBGs could not have been detected. Selection criteria that probe 
large redshift ranges are poorly suited for such overdensity studies. 

We are unable to reproduce the results obtained by K07 who 
found that the 0316 field is a factor 1.5-2 denser in BBGs than 
blank field s. We attribute thi s to th e choice of control field. Ac- 
cording to lvanDokkumetal.1 J2006h the GOODS-S field, used by 
K07, is underdense in 2 < z < 3 red galaxies by a factor of ~2. In 
this study four control fields are used with a total area exceeding 
400n', reducing the effect of cosmic variance on the control fields. 

(iii) The masses and ages of candidate protocluster galaxies 
were determined using SED fitting. LAEs generally have faint con- 
tinuum emission at 3.6 fim and 4.5 pm and relatively blue UV 
slopes, indicating that they contain little dust. Typical stellar masses 
of the LAEs are ~ 10 8 M Q . Four of the LAEs have larger IRAC 
fluxes and stellar masses in excess of 10 10 M Q . This supports re- 
sults in previous studies that there is also a non-negligible number 
of older evolved LAEs 

(iv) The median mass determined for the protocluster LBGs is 
a factor of two larger than that for the field LBGs. Although this 
is not significant, consi dering the uncertain ties of SED fitting, it 
agrees with the results of lLehmer et ail ((2009), who found that pro- 
tocluster LBGs are more luminous at H band and therefore more 
massive than field galaxies. No significant difference between clus- 
ter and field LBGs is found for SFRs and UV-slopes. 

(v) The most massive and intensely star forming galaxies are 
located primarily near to the radio galaxy indicating that proximity 
to the radio galaxy influences galaxy evolution. No other high mass 
and high SFR regions are found in the 03 1 6 field, indicating that the 
radio galaxy is near the centre of the protocluster. The trend in SFR 
is radically different from what is observed in local galaxy clusters 
where the SFR decreases systematically towards the cluster centre. 



We conclude that the protocluster surrounding MRC 03 16-257 
is in a relatively early stage of formation. The blue LBG population 
will likely evolve into a population of passive, red galaxies over a 
timescale of 1 Gyr, forming a structure similar to the well-studied 
protocluster surrounding MRC 1138-262 at z = 2.2 

Spectroscopic follow-up of the LBG and BBG candidates is 
needed to refine their redshifts and establish which individual ob- 
jects are indeed located within the protocluster. Furthermore, popu- 
lation studies of a large sample of radio-selected protoclusters over 
a range of redshift are needed to obtain a more comprehensive view 
of how these structures form and evolve. Such studies would clarify 
whether the protocluster around MRC 0316-257 is representative 
of large-scale overdensities in the early Universe. 
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